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The effective transport coefficients and figure of merit ZT for anisotropic systems are derived 
from a macroscopic formalism. The full tensorial structure of the transport coefficients and the effect 
of the sample boundaries are included. Induced transverse fields develop which can be larger than 
the applied fields and which reduce the effective transport coefficients. A microscopic model rele- 
vant for multi-valleyed materials is introduced which utilizes the effective-mass and relaxation-time 
approximations. The thermopower and Lorentz number are independent of the tensorial struc- 
ture of the transport coefficients in this case and are therefore isotropic. ZT is also isotropic for 
vanishing lattice thermal conductivity K£. For non- vanishing but sufficiently isotropic Ki, ZT is 
maximal along the direction of highest electrical conductivity a. Numerical calculations suggest 
that maximal ZT generally occurs along the principal direction with the largest a j ki. An explicit 
bound on ZT is derived. Several results for specific systems are obtained: (f) Bulk n-type Bi2Tea 
exhibits easily observable induced transverse fields and anisotropic ZTs. (2) Increased anisotropy 
in HgTe/Hgi_ x Cd x Te superlattices (SLs) is associated with larger induced fields. (3) The valley 
degeneracy is split and the bulk masses modified in isolated Bi2Te3 quantum wells, resulting in 
optimal ZTs for wells grown along the trigonal direction. (4) Non-parabolic dispersion in SLs has 
little effect on the thermopower at the carrier concentrations which maximize ZT. 
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PACS numbers: 72.20.Pa, 72.I0.Bg, 73.20.Dx 



I. INTRODUCTION 



In the search to find systems with large thermoelectric figures of merit ZT, the emphasis has been more on new 
materials than-iOa materials structures rQf crystallographic anisotropy. Typical of new structures are superlattices, 
quantum wellsEJu and quantum wires. It has been generally assumed that the direction of highest conductivity in 
an anisotropic material yields optimal thermoelectric properties. The correctness of this assumption is not obvious, 
since directions may exist along which the lattice thermal conductivity is abnormally low and the thermopower is 
high enough to result in large ZT even though the electrical conductivity is less than maximum. 
\ This paper develops a more general transport theory of thermoelectricity in anisotropic systems. The "highest 
■ conductivity" assumption is found to be correct for materials having simple band structures typically of the parabolic 
variety and essentially isotropic lattice thermal conductivities. However, the formalism developed here suggests more 
generally that the optimal orientation corresponds to the principal direction along which the ratio of the electronic to 
lattice thermal conductivity <j/k( is maximum. There are also several surprising results. These include the formation 
of possibly large induced transverse electric fields and temperature gradients, the fact that ZT is strictly isotropic in 
anisotropic systems having parabolic bands if the lattice thermal conductivity is neglected, and that a nearly isotropic 
thermopower and Lorentz number results under these conditions even if the bands are extremely anisotropic and 
non-parabolic. 

The macroscopic formalism, based on the tensorial form of the usual phenomenological transport equations, is 
presented in Sec. |l[ The effects of sample boundaries are included by requiring that the transverse electric and 
heat currents vanish at the transverse surfaces. Anisotropic and isotropic systems are shown to differ both qualita- 
tively, through the presence of induced transverse fields, and quantitatively, through the magnitude of the transport 
coefficients. 



More detailed statements concerning optimal orientations require use of a microscopic model. Section III introduces 
a model commonly used to study transport in semiconductors having multi-valley, anisotropic band structures. The 
transport properties follow from the linearized Boltzmann equation in the effective-mass and relaxation-time approx- 
imations. Intervalley scattering is neglected. The thermopower and Lorentz number somewhat surprisingly turn out 



to be isotropic. The same is true for ZT when the lattice thermal conductivity ki is neglected. Sec. IIIB shows that 



similar results hold for two- and one-dimensional systems. The realistic case corresponding to finite ki is considered 



in Sec. Ill C. If ni is sufficiently isotropic, the maximal ZT is shown to occur for samples cut along the direction 
of highest electrical conductivity. This rather unsurprising result, however, may be modified if the anisotropy of K(. 
exceeds that of the electrical conductivity, leading to the conjecture concerning a j 'n( mentioned above. An explicit 



I 



to 



expression for the upper bound on ZT is also derived, which is a generalization of that previously obtaine' 
anisotropic systems. 

This microscopic description is then applied to several materials of current interest. Sec. IV A examines bulk n-type 
Bi2Te3, in particular, the magnitude of the induced fields, their influence on the effective electrical condu ctivity, and 
the dependence of ZT on sample orientation. HgTe/HgCdTe superlattices (SLs), considered in Sec. IV B, provide an 
example of systems which have a t unable anisotropy. For sufficiently large anisotropy, the induced fields can be much 
larger than the external one. Sec. IV C focuses on isolated n-typc Bi2Tc3 quantum wells. Quantum confinement is 
shown to split the valley degeneracy and modify the effective masses relative to their bulk values. These effects are 
seen to be important in determining ZT quantitatively. 

The effects of non-parabolicity in superlattices on the thermopower and Lorentz number are investigated in Sec. 
Both quantities are seen to be nearly isotropic. The thermopower is also bounded by its values along the principal 
axes, supporting the conjecture that optimal ZT is obtained along those directions. 



II. ELECTRONIC TRANSPORT THEORY 



Consider a rectangular parallelepiped sample of an anisotropic material which is placed between a hot contact at 
x = 0, temperature T/j, and a cold contact at x — L x of temperature T c . An external electric field £ cxt = (£|| Xt ,0, 0) 
is applied to drive a current from the hotter to the colder electrode. Because of the anisotropy, the resulting electric 
current need not be parallel to the x-direction initially; instead it may have a transverse component that leads to an 
accumulation of charge on the faces of the sample, resulting in an induced polarization field. Similarly, the initial 
thermal current through the sample need not be parallel to x, and this may lead to an induced transverse temperature 
gradient. In the following, we compute the effect of these induced fields and temperature gradients on the electrical 
conductivity, thermopower, thermal conductivity, and ZT using phenomenological linear-response theory. We assume 
that the transverse electrical and thermal currents vanish. 

The electric current is 



J e = a{£ - SVT) ee ctF (1) 

where a and S_ are the conductivity and Seebeck tensors respectively. Write F = (F x , F y , F z ) = (Fu , F±) with F» the 
thermoelectric field along x and Fj_ the induced transverse fields. The conductivity tensor takes the form 



T 

QLr,rl & 



(2) 



where crn is the component of the conductivity along the x-direction, and a ± is a 2 x 2 tensor for the transverse 
y, z-directions. Since the full conductivity tensor is not block-diagonal, there will be a 2 x 1 off-diagonal term g_ od and 
its transpose. 

Imposing the boundary condition J e = (J|,0,0) on Eq. (|l|) corresponding to vanishing transverse current gives 

$)-«(&)■ (3) 

Inserting Eq. (|2|) into this expression leads to a linear equation for Fj_ with solution 

F± = -^- 1 <L od F t (4) 
The off-diagonal elements of a induce a current 

T e,ind T xri T —1 771 /r\ 

J\\ = £od F ± = -°-od°-± °-od F \\ (5) 

along the x-direction. Since g_ ± and cr^ 1 have positive eigenvalues, the quadratic form <zJ d g_ ± ~ 1 g_ od is positive definite. 

C.CXt 771 
|| =(T \\ F \\- 



Thus, the induced current opposes the external current J,°' cxt = crui^i. The total current 



II ""II ^"11 

with 



re TC,cxt . T c,ind 

J \\= J \\ + J \\ = fJ l- F ll ( 6 ) 
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The induced fields therefore lead to a reduced conductivity: o\ < a\\. If S_ is diagonal, a\ is the effective conductivity; 
otherwise, the induced temperature gradients give rise to additional terms considered below. In the isotropic case, 
the induced fields vanish, and a\ = an . Note that in general cr 1 > because 

o-iFf = J|fF|| = F T • J e = F T • a ■ F > 0. (9) 

The heat current 3® is given by 

3 Q = TS3 e - kVT, (10) 

where k is the thermal conductivity and the temperature gradient VT = (V^T, V y T, V ' Z T) = (VnT, Vj_T). Using the 
same decomposition as in Eq. (||) for k and S_ in Eq. ( |Io| ) and setting the transverse components of the heat current 
to zero yields 



J?\_f TS\\Ji - K|V T-hJ„,Y_T 



\ I -"II- 1| "II in ^oa 

J ^T^JI-^ViiT-Kj.VxT^ ■ 
Solving for the induced temperature gradients yields 

V ± T = Kl 1 (TS od Jf l -K od X7 ll T). (12) 

Substituting back into Eq. ([ll]) gives 

Jj 3 = TSeff J|f - KeffV||T (13) 

with 

K eS = K|| - K^jftJ 1 /^ (14) 



Observe that oi, Eq. (|7j), and «; e ff; El- (Ml); have the same form. As in the o\ case, the induced temperature 
gradients produce an induced heat current along x which opposes the external heat current — rti|V||T. The net result 
is a reduced effective thermal conductivity: < K c ff < 
Returning to Eq. (|^), we have 

Jf = a^ = G\ (f„ - S||VyT- £ d V ± T) . (16) 
Substituting Eq. @ into Eq. @ leads to 

Jf = <7 eff (£ || - 5offVT||). (17) 

with 

C e ff = y — (18) 

1 + <JiTS_ od K ± S_ od 

and 

Sefi = 5|| - S^kJ 1 ^. (19) 

From the properties of or, < cr e g < mi. Therefore, induced fields and temperature gradients always reduce the 
effective conductivity. Also, S e g in Eq. (|l^) is the same as S e g in Eq. (|l5|) because, by the Onsager relations, k ± and 
hence kJ 1 are symmetric. 

In the steady state, the energy-conservation equation reads 

V-J Q =£-J e . (20) 



3 



Using Eqs. @ and (Q), and assuming that the Thompson heat is negligible, 

V • J Q = V • (TS3 e - kVT) = VT5J e - V • (kVT) (21) 

and 

s ■ r = 3 e { s r 1 3 e + swt) = j e qr x J e + rsyr. (22) 

Thus, Eq. (po|) becomes 

J e £- X J e + V • (kVT) = 0. (23) 

Since J e = (J,?, 0,0), the first term is (ct _1 )iiJ| 2 = (det g_ ± / det cr) J| 2 = erf 1 ^ 2 . The last step is justified by 
expanding the determinants and comparing to Eq. (||). With VT = (V||T, Vj_T) and using Eq. the second term 
is 

= Kcff v 2 r + ^: 1 ^ od rj | f. (25) 



a^jf + « eff VfT = 0. (26) 



Therefore Eq. fl2Q) reduces to 



With the help of Eqs. (0), (13) and (p6|), the thermoelectric figure of merit becomes 

ZT = Ta cS S 2 cS /K cS . (27) 

The transport coefficients are simply replaced by their effective versions. 

These results are derived more succinctly using a general formalism based directly on the Onsager coefficients in 
Appendix A. Rigorous bounds can be placed on the magnitude and sign of <7 e ff and re e ff, but not on S e s without 
introducing a microscopic model that relates these transport coefficients. The next Section considers an example of 
such a model. 

III. MICROSCOPIC MODEL: PARABOLIC BANDS 

A. Three-Dimensional Structures 

According to semiclassical transport theory, the Boltzmann equation in the relaxation-time approximation yields 
the transport coefficients 



a l3 =e z J de{-df /de)J^ (e) (28) 

T{p ■ S)ij = e / tfe(-a/ /ae)Ey(e)(e - /x) (29) 

r«o,« = J de(-fl/ /&)E«(e)(e - rf, (30) 

where /o is the Fermi-Dirac distribution fo(e) — l/(cxp((e — /x)/fcsT) + 1), fi the chemical potential, and 

— Vl (k)v 3 (k)T(k)5(e-e(k)) (31) 

are the components of the transport distribution tensor, the generalization of the function discussed by Mahan and 
SofoO Here e(k) is the electronic dispersion relation, Uj(k) = h^ 1 de(k)/dk i the electronic group velocity, and r(k) 
the relaxation time. Note that k q is the electronic thermal conductivity at zero electrochemical potential gradient 
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inside the sample;0 the usual electronic thermal conductivity at zero electric current, re e , is given in terms of k by 
n e — k — TaS 1 . 

The microscopic model to be used here assumes the conduction to be taking place in a single parabolic band having 
N degenerate valleys centered at k^ n \ n — 1, . . . , N respectively. The dispersion relation for each valley is 

£ («)(k) = £o + (h 2 /2) - fcWjAfW" 1 ^ - kf) (32) 

i,j 

with the inverse effective-mass tensor. The corresponding group velocity is 

v^ik) = h-'de^/dh = f^M^-'ik, - kf ] ). (33) 

i 

Intervalley scattering will be neglected. Thus the transport distribution tensor involves just a sum over the N 
valleys. Assuming the relaxation time to be a function of energy alone, r(k) = r(e(k)), and independent of crystal 
orientation, 

2d 3 k ^ x-^ . ,rwi-i , . ,r^-i . „, . ^u, 



S ij ( £ ) = r(e)X: / f^ 2 E M^-%M^)- 1 k j ,S(e - e(k + k^)) 

n=l i'i' ^ 



3 (35) 



( 27r )3 Z^Z^ JU «' Ju jy 

n— 1 i'j' 
2 3/2 r(£)£ 3/2 

37T 2 n 3 



Q(e-kXS n) k)d 3 k 



(36) 



A" 



J2ydetM {n) Ml ?>-\ (37) 



where are the components of = (fi 2 /2)M^ -1 . This transformation relies explicitly on the validity of the 
effective-mass approximation, the neglect of intervalley scattering, and the relaxation-time approximation used here. 
The square bracket in Eq. ([36]) is evaluated using the identity 

detX (n) = (detX^rf- 1 (38) 



(n) 



and the change of variable k' = \J X^k, Thus, 

E(e) = AT(e) (39) 



with 



N 

A=J2 (mo 1/2 yJdetM (n) ) M^ 1 (40) 



and 



03/2 1/2 

T(e) = » £ 3 / 2 r( £ ). (41) 

The constant, dimcnsionless matrix A contains the full tensorial structure and separates it from the energy dependence 
inT(e). 

The conductivity then becomes 

a = e 2 J de(-df /de)T{e)A = a A. (42) 

Since 
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(vS) = (e/T) J de(-df /de)T(e)(e ~ fi)A = Aa S , (43) 

S = (7- 1 (a-S) =5nl. (44) 
The Seebeck tensor is therefore necessarily isotropic. Further, 

k = Ak (45) 

with 

«o = T- 1 J de(-df /de)T(e)(e - fi) 2 . (46) 

The zero-current electronic thermal conductivity becomes 

«e = - T °±L = A^e (47) 

for n e = kq — To-qSq. 

These results lead to the following surprising conclusion: When the lattice thermal conductivity is neglected, then, 
within the effective-mass approximation as specified here, the thermoelectric figure of merit ZT is independent of the 
sample orientation. Note that 

k e ■ oT 1 = Ak c ■ A^/ao = (K e /<7o)l = L T1, (48) 
where the Lorentz number Lq = n e /aoT. Thus, K e ^j = LfyTo~ij and 

Kcff = L Q To- cB (49) 
as expected. Finally, since S_ is isotropic, S c g = Sq. Combining these effective transport coefficients yields 

ZT = Ta cB S 2 s / k cS = Sq/L , (50) 

a constant independent of direction. 

B. Lower-dimensional Structures 

Dimensionality enters the transport coefficients through the k-space integrals 

d 3 k <-> (2n/L z )d 2 k <-> ({2ir) 2 /L y L z )dk (51) 

for three, two or one dimensions respectively, where L y and L z are the sample sizes in the y- and z-directions. 
Dimensionality also enters through the confinement energies and effective masses for carriers constrained to move in 
a lower-dimensional device. 

For the two-dimensional case with £ — (£° xt , £' nd ) and VT = (V X T, V y T), the induced fields are given by Eqs. (^) 
and (fl2|), the transport coefficients by Eqs. (|l4]), (fl5|), and (|l8| ) with er^ replaced by cr yy , a od replaced by a yx ; similarly 
for the other transport coefficients. The analogue of Eq. (|3l|) for the components of the transport distribution tensor 
is obtained within the effective-mass approximation using Eq. (|5l|): 

Etf (e) = ^rf £ S M®- X M&- X I #U*krS{e ~ ^\). (52) 

z n—l i'j' 

This expression has the same form as Eq. (|3^) with 

N 



A = J2^JdetM {n) Mj n} - 1 (53) 



n=l 

and 



G 



T(e) = eT(e)/nh 2 L z . 



(54) 



These equations are to be contrasted with Eqs. ( |40|) and ( [4l|) for the three-dimensional case. Thus, analogously to 
Eqs. (ji^), (|4^), and (|47]), £ = ^4cro, S_ — SqI and K e = An e where the two-dimensional A and T(e) must be used in 
defining er , So an d K o- Finally, K e cr _1 = L T1 so that 

Kcff = L Ta cS . (55) 

Also, S'cff = So because S_ is isotropic in two dimensions as well. The corresponding figure of merit 

ZT = Ta cS S 2 cS /K cS = S$/L (56) 

is again independent of direction, as in three dimensions. 

It is seen that the two-dimensional and three-dimensional results are entirely analogous and that the former are 
obtained from the latter by taking the limit as one of the effective masses tends to infinity. In this limit, the ellipsoidal 
surface in k-space of constant energy becomes increasingly prolate, until it reaches the edge of the Brillouin zone, 
after which it assumes a cylindrical shape extending from —ixjL z to it / L z upon further increase in the effective-mass 
parameter. Furthermore ZT in two dimensions can be enhanced due to the factor of l/L z in the density of states, 
which, as pointed out in Ref. becomes large for small thicknesses. 

In the one-dimensional case there is no transport in the transverse direction. Thus there are no transverse fields 
and the microscopic and effective transport coefficients are the same. Moreover, all transport coefficients are scalars. 
The transport distribution function is found to be 

N 9 I 9r 

^-ShkVw^ (57) 

Just as in two dimensions, the l/L y L z factor leads to an enhancement of the density of states and thus of ZT for thin 
wires.13 (The wire becomes approximately one-dimensional when it is thin enough that the confinement energy places 
all states associated with its three-dimensional structure sufficiently high in energy that they do not contribute.) 



C. Implications for ZT 

We shall now derive an upper bound for ZT of the Sofo and Mahan formS and show that the highest conductivity 
direction gives optimal values for ZT. The inclusion of an isotropic lattice thermal conductivity Kt causes ZT to lose 
its isotropy. 

The figure of merit including an isotropic lattice thermal conductivity K f — Kgl may be written in the form 

ZT = Ta cS S 2 s / k cS (58) 
where a e g and S e g are as in Sec.g ^eff = k*+ K£ with 

^e^yxi^^e^zz &£) ^e,yz^e,zx 



e,xx ^e.xy , _i_ \( \ \ 

[Ke.yy \ ^i)\^e,zz T" f^i) ~ ^e : yz^e,zy 

^ f^e^zx^^e^yy H~ R>£) ^e,zy^e,yx 

{^e,yy H - K>t){ft>e,zz H~ fc£) ^e,yz^e,zy 



(59) 



defining the electronic thermal conductivity in the presence of the non- vanishing Kg and the sample boundaries. As 
shown in Appendix B, the upper bound on ZT is given by 

ZT < a (K /Ki). (60) 

The dimensionless quantity clq is defined by Eq. (B2). In the isotropic case0 ao = 1; in the present case, clq is of 



order unity. The maximum value of £ defined in Ref. ^ and in Appendix B [Eq. (B12)] of unity is attained if and 



only if T(e) is proportional to a <5-function. In the more physical case Tie) cx e r with r varying between -0.5 and 
2, numerical computations show that £ tends to 1 as /i/kgT — > — oo and to zero as fi/ksT — > oo. For fi/kgT = 0, 
£ ranges from 0.5 to 0.8. Thus, the upper bound can be reached at the cost of going to low carrier concentrations, 
whereas higher carrier concentrations imply smaller £. The optimum concentration lies somewhere in between. 
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We now show that, in the effective-mass, relaxation-time, no intervalley scattering, and isotropic-thermal- 
conductivity approximations, ZT is highest in the direction of maximum electrical conductivity. In the anisotropic 
case, <7 e ff and n c & have an angular dependence obtained in the frame of the sample by a rotation from their common 
principal frame. S c f{ does not because in our microscopic model S_ is isotropic. ZT is therefore also anisotropic. Let 
now P be any symmetric and positive matrix and A > a positive number; then by the properties of positive matrices, 
P^ 1 > (P_ + XI) 1 ■ For the application to thermoelectrics let 



K || Hod 

Hod Hx 



(61) 



In analogy with Eq. (|TJ) we find 

k* = K|| - K od T (K ± + neir 1 ^- (62) 
Taking P = Kj_ and A = «£, the relation P^ 1 > (P + Al) -1 shows that 

Hod T (K±_ + ^ir^od < Hoc? HiT 1 !Hod- (63) 

Thus, for any Ke 

K*e(Ke) > K* e (K( = 0) = (K e /cr Q )(T c ff. (64) 

This shows that the ratio K*/cr e fi is minimized when the axis of current flow in steady state lies along one of the 
principal directions, where, whatever the value of ki, the induced-field related terms vanish and so equality obtains 
in Eq. (jsi]). Writing ZT in the form 

ZT= ; ^ ; , (65) 

K*/Tcr e ff + ne/Ta c s 

the second term in the denominator is seen to be smallest along the principal direction with largest electrical conduc- 
tivity and therefore ZT is maximized for current flow along this direction. 

On the other hand, for a sufficiently anisotropic lattice thermal conductivity the favored direction might be de- 
termined by its minimum rather than that of the highest electrical conductivity. For crystals in which two of the 
principal values of the electrical conductivity are equal, such as Bi2Te3 and SLs, numerical results show that as long 
as the anisotropy of is smaller than that of a, the optimum ZT is still to be found in the direction of greatest 
electrical conductivity. This suggests that, generally speaking, the figure of merit will be maximized in the principal 
crystal direction in which the ratio <j%/kh is greatest, where <7j and Ke t i are the principal values of the electrical and 
lattice thermal conductivity tensors obtained from summing over valleys. 



IV. APPLICATIONS 



In the Section, the formalism of Sec. Ill is applied to several systems in which the effective-mass approximation is 



valid and for which, as shown above, the thermopower is isotropic. 



A. Bulk Bi 2 Te 3 



Bismuth telluride is a semiconductoE-with a room-temperature energy gap of 0.13 eV that crystallizes in a trigonal 
structure with space group Dg d (R3m).c3 An orthogonal coordinate system can be defined in this structure, consisting 
of the trigonal (three-fold rotation) axis, a bisectrix axis that resides within a reflection plane and is normal to the 
trigonal axis, and a binary axis which is along the two-fold rotation axis orthogonal to the other two directions. This 
coordinate system will be used throughout the paper and will be referred to as the crystal frame. 

As shown in Fig. |l|, the constant-energy surfaces of the conduction band at low-doping are ellipsoidal. There are six 
degenerate valleys, each described by a highly anisotropic effective-mass tensor. E3 As illustrated in the figure, two of 
these valleys are bisected by the trigonal-bisectrix plane and their principal axes are rotated approximately 1.2° from 
the bisectrix axis. The light effective masses are mi = 0.025mo near the bisectrix, TO2 = 0.26too along the binary axis 
and 7713 = 0.19mo near the trigonal axis. These two surfaces are related to each other by inversion through the origin. 
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The remaining four constant-energy surfaces are obtained from the first two by rotations through ±27r/3 about the 
trigonal axis. 

The tensorial transport distribution function, Eq. (|39|), involves the scalar part T(e) given by Eq. (|4T|). The 
relaxation time r(e) takep-.to be independent of energy and is determined from the experimental mobility along the 
bisectrix, 1200 cm 2 / V s.liia The matrix pacL^ is constructed from the experimentally derived inverse mass matrix for 
a single conduction band ellipsoid (n = I) S3. Application of the point group operations of the crystal give the other 
inverse mass matrices needed to construct A in the crystal frame [cf. Eq. (f40|)], which is then transformed into the 
frame of the sample. ZT is obtained from the resulting transport coefficients by taking the lattice thermal conductivity 
at 300 K to be 1.5 W / m K (Ref. 13) and isotropic. We shall assume a carrier density n = 5.2 x 10 18 cm -3 , which 
gives the maximal ZT — 0.71 in our model when the sample is cut along the bisectrix. 

The magnitude of the induced fields in a crystal cut at an angle 9 with respect to the trigonal axis in the trigonal- 
bisectrix plane is plotted in Fig. ||(a). The induced field vanishes when the sample is cut along either the bisectrix 
or the trigonal directions. This field reaches its maximum of 76% of the external field for a sample cut 0.367T radians 
with respect to the trigonal axis. As discussed in Sec. y, an external temperature gradient may induce temperature 
gradients along the transverse faces. The induced gradients are at most 12% of the external gradient and have little 
effect on ZT. 

The effective electrical conductivity a e g [Eq. ([l8|)] is shown in Fig. ||b). As the direction in which the sample is cut 
is changed from the bisectrix (9 — n/2) to the trigonal (9 = 0), the conductivities computed with (a c g) and without 
(a xx ) sample boundary effects are both reduced by a factor of four due to the intrinsic anisotropy of the material. 
CT e fi is further reduced with respect to <r xx by the induced fields when the sample is oriented along low-symmetry 
directions by as much as 60%. The influence of the sample boundaries is therefore substantial. 

The combination of intrinsic anisotropy and the effects of the induced fields also affect ZT, but o nly if ki is non- 
zero. As shown in Fig. ||(a) for a hypothetical material with ki = 0, ZT = 2.6 and is isotropic [Eq. (50)]. As seen in 
Fig. [5](b), when K( assumes its bulk value, ZT becomes anisotropic and decreases to a maximum value of 0.71. The 
maxim al ZT applies to samples cut along the high-conductivity bisectrix-binary plane, consistent with the results 
of Sec. [[II C . The minimal ZT = 0.22 occurs for samples cut along the low-conductivity trigonal axis. Despite the 



complicated many- valley band structure, ZT is independent of the azimuthal angle 4> defined in the inset to Fig. [j](a). 
This is consistent with a group-theoretical analysis and is suggested by the relative orientations of the ellipsoids in 
Fig. | 

Introducing anisotropy into kj> by increasing its value along the bisectrix reduces ZT for samples cut along this 
direction, but not for those cut along the trigonal axis. When the ratio a/ ki along the trigonal axis exceeds that 
along the bisectrix, the maximal ZT shifts to the low-conductivity trigonal direction. ZT is ne ver maximal along a 
low-symmetry direction. These numerical results support the conjecture at the end of Sec. Ill C . 



B. HgTe/HgCdTe Superlattices 

We now consider HgTe/HgCdTe SLs, one of many SLs whose well and barrier materials have direct T-point band 
gaps. The constant energy surfaces of the lowest conduction subband CI consist of single ellipsoids of revolution 
aligned with the growth axis of the SL. The conductivity will typically be a minimum along the growth axis (<T m i n ) 
and a maximum within the SL plane (<r max ). If the sample is cut at an angle (j) relative to the SL planes, the induced 
fields will be non-zero only for the faces intersected by the SL growth axis, as indicated in the inset to Fig. ||. Use of 
Eq. (||) then yields the relative magnitude of the induced electric field as 

£ md (c m ax - c min ) sin (f>COS (j) 

T^t = rrr 5— ■ (66) 

CTmaxSm (f> + (T roill COS* (j) 

This relation is plotted in Fig. ^ for several values of the ratio cr min /<7 max . When the system is isotropic, o _ m in/ 'max — 
1 and the induced field vanishes. With increasing anisotropy, £ lnd becomes non-zero away from the high-symmetry 
growth and in-plane orientations and peaks at an angle 

</>max = Sin" 1 y/ Cmin/ (^max + CT min ) ■ (67) 

The maximum induced field is 

£ind^pext| _ ( CTmax _ C r m in)/2 V / (68) 



In the limit <7 m i n /er ma x << 1, these relations reduce to (f> max ~ \/cr m i„/cr m ax and £ md /F cxt | max ~ (l/2)y/a max /cr ri 
Thus, as illustrated in Fig. ||, a relative induced field much larger than unity is obtained for sufficiently small <7 m i n /er m 







From Eq. (|68|), the maximal induced field exceeds the external one when £r m i n /er m ax < 3 — 2\/2 = 0.172. This 
geometrical effect may be practically useful. 

As specific examples, consider two superlattices: (1) 30-A HgTe / 20-A Hg isCdo 85 Te and (2) 20-A HgTe / 40- 
Ah go.15Cdo.85Te. Realistic band structures for these materials are obtained through an 8-band K • p theory within 
the envelope function approximationE-3 based on the parameters in Tab. |. Application of this approach to other 
II- VI and III-V SLs yields-accurate band structures which reproduce experimental optical absorption spectra with 
no adjustable parameters .til The masses of the CI subband along the growth and in-plane directions required for the 
transport calculations are given in Tab. || The room-temperature energy gaps are seen to satisfy E g > lOfcsT in both 
structures, so hole conduction is negligible. Transport in the next highest conduction subband is also negligible, since 
it lies more than 20fcgT higher in energy than CI in both SLs. The relaxation time and lattice thermal conductivity 
enter ZT significantly only in the ratio Ki/cr c sT [Eq. (j65|)]. Their values are estimated from those of an alloy of the 
same composition as the SL (Tab. |J). Thus, interface scattering is neglected. 

The thermoelectric performance of these SLs is shown in Tab. [n| Since a m i n /a max = m^/mj_, SL (l)'s 0.66 mass 
ratio between the in-plane and growth directions yields a maximum induced field of only 0.2F ext . SL (2)'s much 
smaller 0.39 mass ratio, indicating a larger anisotropy, yields a correspondingly larger induced field of 0.5-F . For 
even larger anisotropics, the effective mass approximation breaks down, but estimates using the band structure of 
Sec. suggest that induced fields up to 17F oxt can be obtained in a 40-A HgTe / 50-A CdTe SL. As for bulk Bi 2 Te 3 , 
induced temperature gradients are small, at most 2% of the external gradient for SL (1) and SL (2). Note that the 
maximal ZTs shown in Tab. H although small in magnitude, represent a substantial increase over the equivalent bulk 
alloy. 



C. Bi 2 Te 3 Quantum Wells 

Consider a quantum well formed from a multi- valley semiconductor. Within the single-band envelope function 
approximation,!!^ the wave function for the nth valley, ip^ (r), satisfies 



-5>g l) ^+y(r) 



V> (n) (r) = £^V (n) (r), 



(69) 



where a l} = (h 2 /2)M^ ) ' is proportional to the bulk inverse effective mass tensor M^f \ E^ the energy measured 
from the bottom of the valley in the bulk, and i, j — x, y, or z. V(r) is the confining potential, which is taken to be 
a square well of width L z and infinite height. Solution of Eq. fl69|) yields the eigenfunctions 



( r ) = V 2 / L z exp i(k x x + k y y- (a^k x + a^k v )z/a(™A sm(mirz/L z ) 



(70) 



and the eigenenergies 



with 



, (n)Well, 



(71) 



(n)Woll 



M 

CXxx 
(«) 



(n) (n) ; In) 
Oixz a Z x J Otzz 
(n) 

ZZ 



Oizx jOi 



(n) 

a xy 



/„,(") 

<-Xxz 

ZZ 



tXzy j tXz 
j a 



(72) 



and the confinement energy 



^(n)Co 



a 



Wn 2 m 2 /L 2 z , 



(73) 



for m = 1, 2, 3, ... 

Thus, the quantum well band structure differs from that of the bulk in two ways. First, the energy of the bottom 
of the rtth valley is shifted with respect to the bulk by an amount £;( n ) Conf that is inversely proportional to m^z- 

(n) 

Since m\z is generally different for different n, the confinement splits the valley degeneracy. If the well is sufficiently 
narrow, only the valleys having the lowest confinement energy will contribute to the transport. Second, as seen from 
Eq. ([72]), the masses in the quantum well differ from those in the bulk by terms second order in the off-diagonal 
elements of the bulk inverse effective-mass tensor. 
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Suppose the well material is n-type Bi2Te3. The bulk band structure was described in Sec. IVA. The relaxation 
time and lattice thermal conductivity can be obtained as in Ref. [|, which neglects interface effects. Eqs. (|72]), (|39|), 
and (|5^)-([34]) then lead directbi-|to the transport coefficients and ZT. For simplicity, m of the barriers is neglected. 
Thus, ZT is an overestimate.ErE£l 

Examine two orientations of the quantum well: one in which the growth direction is along the trigonal axis and 
one in which it is along the binary. In both cases, take the currents to flow along the bisectrix. These choices are 
expected to yield the best ZTs, since the currents are along the bulk high-conductivity direction and the alignment 
of the well with the crystal axes eliminates detrimental induced fields. The maximal ZT for these configurations is 
plotted as a function of well width in Fig. |[ 

As seen in this figure, quantum well growth along the trigonal axis yields larger ZTs. From Fig. [l], the masses of each 
ellipsoid along the trigonal axis are the same, but they differ along the binary axis. Thus, all six ellipsoids contribute 
to transport in wells grown along the trigonal axis, but valley splitting allows only two ellipsoids to contribute in wells 
grown along the binary axis. The resulting larger density of states in the former case leads to higher ZTs. For a given 
orientation, the density of statesjs also larger for smaller well widths, resulting in an increase of ZT with decreasing 
well widths in both orientations .u 

For comparison, Fig. |^ also shows ZT computed using the simpler Bi2Te3 band structure of Ref. ^. This band 
structure consists of six equivalent ellipsoids whose principal axes are all aligned with those of the crystal. As in the 
multi-ellipsoid model, ZT is larger for wells grown along the trigonal axis and for smaller well widths. However, the 
magnitude of ZT is larger than in the multi-ellipsoid model. This difference arises because each ellipsoid shares the 
same orientation, allowing all of the ellipsoids' low-mass directions to lie along the bisectrix and preventing any valley 
splitting. 



V. MICROSCOPIC MODEL: NON-PARABOLIC BANDS 



Within the effective-mass approximation, the analysis of Sec. Ill showed that the thermopower and Lorentz number 
are isotropic. To see how a noB-parabolic band structure affects these conclusions, consider an electronic dispersion 
relation of the Esaki-Tsu formO 



e(k) = 7T fcjj /2m\\ + A(l - cos k z d) (74) 

with wave vector k = (fc|| cost/?, fc|| sin (p, k z ). This relation models a supcrlatticc of period d. The in-plane dispersion 
is parabolic with mass mn, and that along the growth direction has a tight-binding form with band width 2A = 
(e(0,0, 7r/d) — e(0, 0,0)). Since the mass along the growth direction m z — h 2 /Ad 2 , the anisotropy can be increased 
by reducing A. 

In the principal frame of the SL, the transport distribution tensor [Eq. (|3l|)] is diagonal with components 



V ,,,( c) - m W dT ^ { A 2 cos- 1 (l- £ /A) + ( £ -A)Ve(2A- £ ), e < 2A 

^Growth(£)- < nA2 ^ £ > 2A (75) 



along the growth direction and 

s (] _ r(e) f ( £ - A)cos- 1 (l-e/A) + ,/ £ (2A - £ ), £ < 2A , , 



along the planes. The transport coefficients are obtained from Eqs. (p8|)-(p0|) in the principal frame and then trans- 
formed into the sample frame. As in Sec. H, the relaxation time r(k) is assumed to be a function of energy only. 
Direct computation with r oc e r indicates that the qualitative features discussed below are independent of the choice 
of r. Quantitatively, the thermopower is an approximately linear function of r at fixed sample orientation and chem- 
ical potential /x and increases by approximately 50% as r goes from to 1.5. In what follows, r — 0, d — 100 A, 
m||/mo = 0.021, and A = 57 meV, corresponding to the CI subband in a 50-A Hgo.75Cdo.25Te / 50-A Hgo.7Cdo.3Te 
SL. 

The anisotropy in the resulting effective thermopower S c e [Eq. (|l5|)1 and Lorentz number Lo^s = K>efi/&e&T is 
shown in Fig. ||. For /i < 0, S c g along the growth and in-plane directions differ by < 10%. This near isotropy is 
expected, since the carriers determining S c g are near the zone center, where the effective-mass approximation is good. 
The anisotropy increases substantially as \x increases past 2A, reaching over 6000 at fi = 0.4 eV. From Fig. ^(b), 
the anisotropy in Lo, e ff is at most 30% and goes to zero in the largc-/i limit, approaching the metallic value of 
(tt 2 /3)(fce /e) 2 . As ZT for this band structure is maximal for /j, *~ 0, the anisotropy in S e g and Lo.eff in the relevant 
parameter range are small. 
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These numerical results are also consistent with the idea that the induced fields are detrimental to ZT, as discussed 
in Appendix A. As shown by the dotted lines in Fig. |^(a), the induced fields present along low-symmetry directions 
reduce S e s below S xx . However, S c s is always bounded by S c g = S xx along the principal directions. Thus, optimal 
thermoelectric performance is achieved for samples cut along the principal axes. 



VI. SUMMARY 



This paper developed the transport theory relevant for anisotropic, multi-valleyed materials, taking into account the 
full tensorial structure of the electronic transport coefficients and including the effects of sample boundaries. Induced 
transverse fields are associated with samples cut along directions in which the conductivity had off-diagonal elements. 
These fields can be larger than the applied fields. Within the effective-mass and relaxation-time approximations and 
the neglect of intcrvalley scattering, the tensor characterizing the structure appears as a simple multiplicative factor 
of each of the transport coefficients. This factorization results in an isotropic thermopower and Lorentz number even 
in the presence of anisotropy. In a hypothetical material having no lattice thermal conductivity ZT would therefore 
be isotropic. For non-vanishing and sufficiently isotropic lattice thermal conductivity ZT is maximal for samples 
cut along the high- conductivity directions. Widely ranging numerical results suggest that the maximal ZT generally 
occurs along the principal direction which maximizes the ratio of the electronic to lattice thermal conductivity. An 
upper bound for ZT is given which generalizes the result of Mahan and SofoEiJ to anisotropic systems. 

Application of these formal results to some systems of interest showed the following. (1) Bulk n-type Bi2Te3 
exhibits induced transverse fields, reduced effective conductivity, and anisotropic ZTs. These effects should be easily 
observable. (2) In HgTe/Hgi_ x Cd x Te SLs, whose anisotropy is adjustable by varying the layer widths and composition, 
increased anisotropy is associated with larger induced fields. (3) The confining potential in isolated Bi2Te3 quantum 
wells splits the valley degeneracy and changes the masses from their bulk values. The latter features suggest that 
optimal ZTs are associated with wells grown along the trigonal direction. (4) The effect of non-parabolic dispersion 
on the thermopower is very small for the doping range over which ZT is maximal. The maximal thermopower occurs 
along principal directions, which supports the suggestion given in the previous paragraph. 
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APPENDIX A: 



In this appendix we present a unified form of the linear-response theory which, because of its explicit use of the 
Onsager coefficients, has the advantage of treating the electric and heat currents more symmetrically than in Sec. |J. 
Note that, since the Onsager coefficients are assumed to be symmetric matrices, this formalism does not apply to 
monoclinic, triclinic, or chiral materials. The generalized fluxes J are given in terms of the generalized forces f by 



with 



and 



Lf (Al) 



Jx Jx Jy J*y Jz Jx ) — ( T " ) (A2) 



f=(e£ x -W X T/T e£ y ~V y T/T e£ z -V Z T/T f = j (A3) 

where Jy and fy contain the first two components and and fj_ the remaining four components of J and f , respectively. 
With respect to these conventions the matrix of Onsager coefficients, L, assumes the form 



12 



with 



in 



L 



-21 



j 11 r: 
12 r22 



(A5) 



L 



•od 



T 11 

^xy 

12 



7-21 

22 



-21 



rll t: 

xz xz 

t 12 r22 rl2 r 22 

xy xy xz xz 



(A6) 



and 



L ± 



Equation (Al) and the condition that J± ~ yield 







t21 

yy 


T 11 

^yz 


T 21\ 

L yz \ 






j 22 


T 12 


7-22 




yy 


yy 


^yz 








T 21 

^zy 


LH 


Lf z 




\ rl2 


t22 

L zy 




Lfz) 


= yield 










f± = 


-L.1 1 


Liodh ; 





(A7) 



hence, 



with 



J|| - i c ff f || 



I^ S = L {l -L^ d Ll 1 L od ^L ll ~N ll 



It can be shown that N_» > and i cff > 0. We may decompose Jii into external and induced parts as follows: 



ext i -rind 



with 



and 



Jii = J 



J oxt =Liifii 



rind 



= -iViifn 



(A8) 
(A9) 
(A10) 

(All) 
(A12) 
(A13) 



In Sec. H it was found that, for an applied electric field only (f|| = (e£ x ,0)), the induced electric current opposes the 
external electric current because a c g < au , and likewise for an applied thermal gradient only, the induced heat current 
opposes the external heat current. Physical intuition suggests that these results should generalize to the case where 
both an applied electric field and applied thermal gradient are present. Then we expect that, if both components of 
f have the same sign, J md should oppose J cxt in the sense that they lie in opposite quadrants of the plane spanned 
by J| and J§ . Although we have not proven this conjecture, it is consistent with numerical results for randomly 
constructed matrices of Onsager coefficients. Undoubtedly a proof requires a general microscopic model. Finally, we 
note that the entropy generation is given by 



dS 
-T — 
dt 



J f 



•fii 



■fii. 



(A14) 



field 



Since L« > and N_« > 0, we have J cxt • f|| > and J • f|| < 0. Thus, external currents decrease S|fi e id but induced 
currents increase S , |fi c id- 
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APPENDIX B: 

We present here a derivation of an upper bound of the Sofo and Mahan forrnQ for an anisotropic material. From 
a = Axjq we infer 

c c ff = cr a (Bl) 

with 

A A — A A A A — A A 

, t ^yx^zz -"lyz-^zx . -^zx-^yy -^zy-^yx /oo\ 

a = A xx - A xy - - — — — A xz — — . (UZ) 

yy zz y z z y zz yy y z z y 

When the K^-dependence in Eq. (|5^) is taken into account, we find similarly 

k* = n e a(Ki) (B3) 

with a(ftf = 0) = ao and 



— A xx — A xy 



Ay X {A zz + ki/ k c ) Ay Z A zx 

{Ayy + Kl/ K e )(A ZZ + K£/K e ) ~ Ay Z A Z y 



A zx (A yy + Ke/K e ) — A zy A yx 

{Ayy + Kl/ K e )(A ZZ + Kg/ K e ) — Ay Z A Z y 



Therefore 

Ta S 2 a 
(k - Ta S%)a(Ke) + K e 

Following Mahan and Sofo, introduce dimensionless integrals 



^=7-S^-- (B5) 



/ e 
dx / x , ^o s ( x ) xn i s{x) = hr T((i + xkT), (B6) 
-M/fesT (e x + 1) 

where ro is the Bohr radius. In terms of these moments, 

(TqSq = {k B /e)a Q Ii (B7) 
«o = {k B /e) 2 Td h (B8) 

where CTq = e 2 /hr has dimensions of conductivity. Then 

ZT = T(k B /e) 2 a 2 I 2 /a a Io 

((k B /e) 2 Ta l2 ~ T{k B /efd Q Il/Iv)a{K t )/a Q + K e /a V ; 



(a/ 2 - ctlf / I )a(Ki) / a Q + l/a 



(BIO) 
(BH) 



(1 — £)a(«^)/ao + .B 
with a = (k B /e) 2 Tao/ Kg, 

e = AW2, (B12) 

and B = l/a^aa — ki/kqciq. By the Cauchy-Schwarz inequality, < £ < 1. The limit as £ t ends to 1 ma xim izes the 
figure of merit by maximizing the numerator and minimizing the denominator in Eq. (Bll). From Eq. (B4) it may 
be seen that for ki/k c — Kg/ (k B / e^Ta^^iX ~ £) sufficiently large, cl(ki) tends to A xx , and certainly this condition 
is met as 1 — £ tends to zero. Thus as £ — > 1 we have that 

ZT^- ^L- <I = ^. (B13) 

(1 - OAw/ao + B ~ B ki v ' 
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TABLE I. Room temperature energy gap E g , spin-orbit splitting A, and heavy-hole effective mass tuhh for a 
HgTe/Hgo.i5Cdo.85Te superlattice. The valence band offset is 0.30 eV from a linear interpolation of the value in Ref. |is| , 
and the momentum matrix element (2/?no) < Z\p z \S > | 2 = 18.0 eV (Ref. po| ). 



HgTe Hgo. 15 Cdo. 8 5Te 



E g at 300 K (eV) a 


-0.166 


1.191 


A (eV) 6 


1.0 


0.9 




0.7 


0.7 



a Rcf. 18 
b Ref. 21 
c Ref. 22. 



TABLE II. Parameters and thermoelectric performance in HgCd/HgCdTe superlattices (SLs). m± and my are the CI 
effective masses along and normal to the growth axis, (£ md /-F cxt ) max is the maximal relative magnitude of the induced field, 
'/'max is the angle relative to the SL planes at which the sample must be cut to obtain this field, (ZT) max is the maximal figure 
of merit for the SL, and (ZT) a n oy is that for the equivalent alloy. 





SL(1) 


SL (2) 




30-A HgTe / 


20-A HgTe / 




20-A Hgo.15Cdo.85Te 


40-A Hgo.15Cdo.s5Te 


Band structure:" 






E g (eV) 


0.25 


0.46 


mx/mo 


0.038 


0.110 


my/mo 


0.025 


0.043 


m\\/m± 


0.66 


0.39 


Transport parameters: 






r(s) b 


5.3 x 10~ 14 


2.6 x 10~ 14 


Ki (mW / cm K) c 


8.2 


12 


Thermoelectric performance:^ 








0.21 


0.48 


0ma X (rad) 


0.22tt 


0.187T 


(^T) m ax 


0.14 


0.086 


(ZT) a ll y 


0.12 


0.061 



"Based on the parameters in Tab. [I]. 
b From the mobility data presented in Ref. 23 



c Quadratic fit to the data in Ref. 24 
d This work. 
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FIG. 1. Constant-energy surfaces of the conduction band valleys in Bi2Te3. The bisectrix, binary, and trigonal directions 
are indicated. 



FIG. 2. Effect of induced transverse fields in bulk n-type Bi2Te3. (a) Normalized induced field F z /F x [Eq. and (b) 
microscopic (a xx , dot-dashed line) and effective (cr e g, solid line) conductivities [Eq. (]l8|)] for the sample geometry shown in the 
inset. 

FIG. 3. ZT of n-type Bi2Te3 at 300 K for (a) a hypothetical material with m = and (b) ki = 1.5 W / m K, the bulk 
in-plane value, as a function of the direction in which the sample is cut (defined in the inset). 

FIG. 4. Relative magnitude of the induced electric field £ md /^? cxt m a superlattice as a function of the angle <j> defined in 
the inset for several values of the conductivity ratio o- min /cr max [Eq. ([j6|. 



FIG. 5. Maximal ZT at 300 K for an isolated Bi2Te3 quantum well within the multi-ellipsoid model of Sec. IV C as a 
function of well width for wells grown along the trigonal (heavy dot-dashed line) or binary (heavy solid line) directions. Results 
for the equivalent-ellipsoid model of Ref. |l for wells grown along the trigonal (thin dot-dashed line) and binary (thin solid line) 
directions are also shown. 

FIG. 6. (a) Effective thermopower S c s and (b) effective Lorentz number Lo.eff = i^cB/<JesT as a function of chemical 
potential for the superlattice described in Sec. |^. Shown are samples cut along <j> = (heavy solid line), 7r/4 (heavy dotted 
line), and 7r/2 (heavy dashed line) as defined in the inset to Fig. [|. The results for cf> = 7r/4 neglecting the effects of the induced 
fields are also presented (light dotted line). The metallic limit of Lo/(fcs/e) 2 = 7r 2 /3 is indicated in (b). 
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